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ABSTRACT 


An  understanding  of  ocean  acoustic  fields,  their  statistics,  and  relation  to  the 
oceanographic  environment  is  the  sine  qua  non  of  undersea  warfare.  In  the  tactically 
important  Philippine  Sea,  powerful  mesoscale  eddies  can  have  strong  effects  on  acoustic 
fields.  To  quantify  eddy  effects,  a  mesoscale  sound-speed  model  was  developed  and 
interfaced  with  a  parabolic  equation  acoustic  simulation.  Eight  combinations  of 
frequency  (20Hz/250Hz),  wavenumber  spectra  (Stammer/Lorentzian),  and  source  depth 
(50m/200m)  were  simulated  through  the  model.  For  each  combination,  the  unperturbed 
transmission  loss  (TL)  curve  and  composite  eddy-field  TL  curve  were  compared  to  assess 
acoustic  variability  caused  by  mesoscale  ocean  features.  Eddies  alter  acoustic  energy  by 
shifting  convergence  zones,  driving  energy  into  the  seabed,  trapping  energy  in  surface 
ducts,  and  increasing  scintillation.  These  effects  are  greater  at  higher  frequencies  and 
deeper  source  depths,  shifting  both  the  mean  TL  difference  and  RMS  variability  on  order 
5-lOdB.  The  wavenumber  spectrum  showed  no  significant  effect  on  acoustic  variability. 
Eddies  also  cause  horizontal  out-of-plane  scattering.  Ray  equations  were  manipulated  to 
demonstrate  that  eddy-induced  bearing-angle  errors  can  be  on  order  one  degree  at  500km, 
increasing  as  the  square-root  of  range.  Target  localization  errors  due  to  angle  error  are  on 
order  7km  at  500km,  increasing  as  range  to  3/2  power. 
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INTRODUCTION 


The  kinetic  energy  associated  with  ocean  circulation  is  99%  contained  within 
mesoscale  ocean  features  having  horizontal  spatial  scales  less  than  100km  [Munk  et  al. 
1995],  Ocean  mesoscale  processes  include  fronts,  jets,  filaments,  squirts,  eddies,  and 
meanders  [Stammer  1997],  Variability  in  kinetic  energy  can  lead  to  significant 
differences  in  sound  speed  through  seawater.  Ocean  eddies,  one  specific  family  of 
mesoscale  features,  produce  sound  speed  fluctuations  of  order  5-15  m/s  root  mean  square 
(RMS)  that  are  superimposed  on  the  background  sound  speed  profile  [Munk  et  al.  1995]. 
Typical  scales  of  these  eddies  are  hundreds  of  meters  in  the  vertical,  hundreds  of 
kilometers  in  the  horizontal,  and  weeks  to  months  in  duration  [Stammer  1997,  Kobashi 
and  Kawamura  2001].  As  such,  these  features  can  lead  to  important  variations  in  acoustic 
fields  to  be  used  in  signal  detection,  classification,  and  localization. 

A.  PREVIOUS  WORK  AND  LINKS  TO  FUTURE  PROJECTS 

Over  sixty  years  of  ocean  acoustic  research  exists,  prompted  by  the  discovery  of 
the  deep  sound  channel  in  1944  [Worcester  and  Spindel  2005].  Observationally,  eddy 
effects  on  acoustic  variability  were  documented  as  early  as  1973  [Vastano  and  Owens 
1973],  but  present  theory  and  computer  modeling  are  quite  limited  [Smith  et  al.  1992, 
Munketal.  1995]. 

More  recently,  the  North  Pacific  Acoustic  Laboratory  (NPAL)  commenced 
research  and  experimentation  in  the  field  of  acoustic  variability  as  influenced  by  the 
mesoscale  environment.  These  and  other  experiments  have  led  to  some  quantification  of 
eddy  effects  on  propagation  loss  and  horizontal  diffraction  [Kirby  1988,  Weinberg  and 
Clark  1980],  but  a  theoretical  basis  for  understanding  these  effects  is  still  insufficient 
[Worcester  and  Spindel  2005].  Further  work  in  prediction  of  eddy  effects  has  been 
accomplished  by  Comuelle  and  Howe  [1987];  Smith  et  al.  [1992];  and  Munk  et  al. 
[1995]. 

Comuelle  and  Howe  [1987]  studied  eddy-induced  ray  loop  resonances  and 
demonstrated  that  eddies  with  wavelengths  near  the  ray  loop  length  have  the  greatest 
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effect  on  acoustics.  Smith  et  al.[1992]  examined  ray  chaos  as  it  applies  to  eddy  effects. 
This  research  suggests  that  deterministic  predictions  for  near-axial  acoustic  energy  using 
ray  theory  are  valid  at  ranges  less  than  one  to  two  thousand  kilometers.  Munk  et  al. 
[1995]  discussed  ray  loop  resonance,  as  well  as  horizontal  angle  deviation  through  a 
single  theoretical  eddy. 

Up  to  this  point,  previous  work  has  not  had  much  emphasis  on  transmission  loss 
curves  or  three-dimensional  acoustic  processes,  such  as  horizontal  out-of-plane  refraction 
caused  by  ocean  eddies.  This  study  is  designed  as  a  small  part  of  the  larger  upcoming 
field  exercise  by  NPAL  in  the  Philippine  Sea.  Thorough  understanding  of  the  expected 
ocean  variability  due  to  the  mesoscale  environment  will  allow  better  development  of  the 
experimental  plan.  Additionally,  quantification  of  the  eddy-induced  acoustic  intensity 
change  measures  the  value  of  accurately  mapping  the  mesoscale  ocean  environment. 

The  most  current  version  of  the  NPAL  Philippine  Sea  Experiment  plan  as  of  April 
2008  is  designed  around  a  six-node  transceiver  array  for  acoustic  tomography,  as  shown 
in  Figure  1 .  The  experiment  aims  to  characterize  the  barotropic  and  baroclinic  structure 
of  the  ocean  environment,  spatially  and  temporally,  in  order  to  generate  an 
eddy-resolving,  four-dimensional  sound  speed  field  for  acoustic  predictions.  The 
scientific  goals,  from  the  Strawman  Philippine  Sea  Experiment  Plans  [working  document 
from  NPAL  11*  Annual  Workshop,  April  2008],  are  to: 

(i)  understand  the  impacts  of  fronts  and  eddies  on  acoustic  propagation  in 
this  highly  variable  region,  (ii)  determine  whether  acoustic  methods, 
together  with  satellite,  glider  and  other  measurements  and  coupled  with 
ocean  modeling,  can  yield  estimates  of  the  time-evolving  ocean  state 
useful  for  making  improved  acoustic  predictions  and  for  understanding  the 
local  ocean  dynamics,  (iii)  improve  our  understanding  of  the  basic  physics 
of  scattering  by  small-scale  oceanographic  variability  due  to  internal 
waves  and  spice,  and  (iv)  characterize  the  ambient  noise  field,  particularly 
its  variation  over  the  year  and  its  depth  dependence. 
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Figure  1.  Array  node  locations,  as  currently  proposed,  for  upcoming  NPAL  field 
experiment  in  the  Philippine  Sea  in  2009-2011.  Color  bar  denotes  underwater 
bathymetry.  Contours  from  500m  to  5000m  are  shown,  at  500m  intervals. 
[Strawman  Philippine  Sea  Experiment  Plans,  working  document  from  NPAL  11* 
Annual  Workshop,  April  2008] 


B.  MOTIVATION  AND  STATEMENT  OF  RESEARCH 

Thoroughly  understanding  the  ocean  acoustic  fields,  their  statistics,  and  their 
relation  to  the  oceanographic  environment  is  the  sine  qua  non  of  undersea  warfare. 
Eddies  are  dynamic  features  which  can  cause  large  sound  speed  fluctuations  because  they 
have  different  physical  properties  than  the  ocean  around  them.  These  localized  acoustic 
fluctuations  on  order  lOdB  can  have  significant  effects  on  target  detection  and  targeting 
accuracy.  Observationally,  physical  processes  by  which  eddies  affect  transmission  loss 
are:  causing  increased  scintillation,  shifting  convergence  zone  ranges,  driving  energy  into 
the  seabed,  and  forcing  energy  into  a  surface  duct. 
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Eddies  are  unlike  internal  waves  in  that  they  ean  be  modeled  deterministieally. 
Theoretieally,  these  deterministie  models  will  allow  scientists  to  predict  the  magnitude  of 
eddy-induced  sound  speed  variability. 

Additionally,  eddies  refract  sound  out  of  the  horizontal  plane,  which  affects  beam 
forming  capabilities  and  accuracy  of  bearing  angle  measurements.  Previous  research  has 
shown  horizontal  refraction  through  mesoscale  eddies  and  fronts  ranging  from  0.5  to 
several  degrees  [Munk  1980,  Weinberg  and  Clark  1980]. 

In  order  to  further  investigate  the  points  listed  above,  the  objectives  of  this  study 

are  to: 

•  identify  the  important  eddy  scales, 

•  assess  the  feasibility  of  constructing  a  useful  mesoscale  model  based  upon 
present  observations  that  can  be  used  for  acoustic  prediction  calculations, 

•  quantify  the  effects  of  eddies  on  transmission  loss  estimates  based  on 
numerical  simulation, 

•  analyze  the  importance  of  out-of-plane  scattering  effects. 

Chapter  II  describes  basic  principles  used  in  creating  the  acoustic  model.  Chapter 
III  discusses  methods  of  analysis,  including  the  parabolic  equation  method,  the  boundary 
conditions  for  the  computational  domain,  and  modeling  of  the  ocean  background  and 
eddy  fields.  Chapter  IV  addresses  the  results  of  the  computer  simulations,  and  Chapter  V 
provides  some  conclusions  and  statements  regarding  future  work. 
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II.  BASIC  PRINCIPLES 


A.  ROSSBY  WAVES 

In  this  study,  eddies  are  modeled  as  superpositions  of  linear  Rossby  waves,  as  in 
the  treatment  by  Munk  et  al.  [1995],  This  method  is  convenient  in  that  it  allows 
computationally  efficient  modeling  of  non-linear  eddies  using  linear  methods  that 
produce  similar  temporal  and  spatial  scales.  As  shown  in  Figure  2,  Rossby  wave 
dynamics  are  unique  in  that  the  dispersion  relation  is  circular: 

cj=-pkj{K^+k;+x^.  (1) 

Here,  <j  is  the  Rossby  wave  frequency,  /?  is  the  latitudinal  gradient  of  Coriolis  force,  and 
k^  and  k^  are  horizontal  wavenumbers.  The  Rossby  eigenvalue  is  represented  by  X , 

where  ^  is  the  Rossby  radius  of  deformation  (the  length  scale  at  which  rotational 
effects  and  buoyancy  effects  are  equally  important). 

High  frequency  Rossby  waves  have  small  wavenumbers  and  long  wavelengths, 
while  low  frequency  Rossby  waves  have  large  wavenumbers  and  short  wavelengths. 
Low  frequency  (less  than  one  cycle  per  year),  large  wavenumber  Rossby  waves  are  used 
here  to  represent  mesoscale  eddies  that  have  the  most  effect  on  acoustics,  because  they 
are  scaled  near  the  ray  loop  length  [Comuelle  and  Howe  1987]. 
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Figure  2.  Rossby  wave  dispersion  relations  for  <7=0.5,  1,  and  2  cycles  per  year. 
The  shaded  regions  represent  those  Rossby  waves  that  most  affect  acoustic 
propagation. 


B.  GLOBAL  WAVENUMBER  SPECTRA 

Until  Stammer’s  work  in  1997,  eddy  dynamics  had  not  been  well  understood,  due 
in  part  to  the  inadequacy  of  available  ocean  observing  systems.  However,  the  forcing 
mechanisms  for  eddies  and  the  methods  by  which  eddies  disperse  energy  are  particularly 
important  to  understand  because  99%  of  the  ocean  kinetic  and  potential  energy  is  carried 
and  redistributed  through  eddies  [Stammer  1997]. 

Satellite  altimetry  provides  a  global  observation  system,  and  therefore,  can  be  a 
valuable  tool  in  gaining  insight  into  mesoscale  ocean  features.  Stammer  [1997]  analyzed 
three  years  of  sea  surface  height  (SSH)  data  from  the  TOPEX/Poseidon  (T/P)  satellite 
altimeter  and  derived  a  universal  spectral  model  of  SSH  to  infer  scales  of  ocean 
variability  in  specific  ocean  regions. 
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The  linear  relationship  between  the  first-mode  Rossby  radius  of  deformation  and 
the  eddy  scale  suggests  that  first-mode  processes  dominate  SSH  variations.  Stammer 
derived  a  global  spectral  relation  for  SSH,  displayed  in  Figures  3  and  4a: 
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by  normalizing  the  along-track  wavenumber  by  the  local  eddy  wavenumber  of  maximum 
k/ 


energy  ( A:  =  %  ). 
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Figure  3.  Zonal  averaged  wavenumber  spectrum  normalized  by  ko,  the  local 
wavenumber  of  maximum  energy.  (After  Stammer  1997.) 

Stammer’s  results  indicate  baroclinic  instability  is  the  dominant  eddy  source  in 
the  extratropical  oceans.  Within  the  framework  of  baroclinic  instability,  the  Rossby 
radius  of  deformation  is  closely  related  to  the  spatial  scale  of  maximum  eddy  growth 
rates.  In  observation,  non-linear  processes  act  to  increase  the  length  scales  of  the  most 
unstable  waves. 


In  addition  to  the  spectrum  designed  by  Stammer,  another  power  spectrum, 
colloquially  named  the  Lorentzian  spectrum  and  shown  in  Figure  4b,  is  another  method 
to  describe  the  wavenumber  relation  for  the  mesoscale  ocean: 


^  k^  k  if 


2  ’ 


W=' 

where  kif  =  l/300km  and  k^  =  k^,ky  are  horizontal  wavenumbers  [Munk  et  al.  1995]. 


(3) 
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Figure  4.  (a)  Modeled  wavenumber  eddy  spectrum  as  designed  by  Stammer  [1997], 

where  frequency  is  a  power  function  of  wavenumber.  Power  values  for  each 
slope  are  indicated,  (b)  Modeled  Lorentzian  wavenumber  eddy  spectrum 
[Munk  et  al.  1995],  where  frequency  is  a  power  function  of  wavenumber  at  the 
higher  values  of  wavenumber,  but  rolls  off  at  k*  as  indicated. 

C.  VERTICAL  PARTITIONING  OF  KINETIC  ENERGY  IN  THE  OCEAN 

Sea  surface  altimetry  provides  valuable  information  in  quantifying  ocean  energy 
at  the  surface;  however,  understanding  the  vertical  structure  of  the  ocean  is  essential  to 
properly  assessing  the  altimetry-derived  sea  surface  elevation  and  slope.  Wunsch  studied 
two  decades  of  moored  current  meter  observations  from  many  sources  and  devised  a 
reasonable  first  estimate  of  vertical  partitioning  of  the  ocean.  Further  work  is  still 
necessary  due  to  the  disparate  spatial  and  temporal  nature  of  the  data;  however  these 
results  are  a  good  starting  point  [Wunsch  1997]. 

For  the  main  portion  of  the  North  Pacific  Ocean  (bounded  by  135“  and  245“E  and 
25°  and  60°N),  about  35%  of  total  kinetic  energy  is  contained  in  the  barotropic  mode, 
with  the  remaining  65%  in  baroclinic  modes.  In  the  area  of  the  NPAL  experiment, 
approximately  30%  kinetic  energy  is  contained  in  the  barotropic  mode,  40%  in  the  first 
baroclinic  mode,  10%  in  the  second  baroclinic  mode,  and  the  remainder  in  higher  modes. 
In  terms  of  surface  kinetic  energy,  as  much  as  90%  of  the  energy  is  contributed  from  the 
baroclinic  modes  in  parts  of  the  North  Pacific. 
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These  findings  are  based  on  searee  data,  few  degrees  of  freedom,  and 
considerable  subjectivity  in  contouring;  however,  they  are  validated  by  altimeter 
measurements.  The  extrapolations  of  current  meter  modes  to  surface  kinetic  energy 
match  up  fairly  well  compared  to  the  altimetric  surface  kinetic  energy  distribution.  A 
useful  outcome  of  the  NPAL  experiment  will  be  a  quantitative  measure  of  proportional 
mode  energy  in  Philippine  Sea  eddies. 
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III.  METHODS  OF  ANALYSIS 


A.  PARABOLIC  EQUATION  METHOD 
1.  The  Parabolic  Equation 

Introduced  to  underwater  acoustics  in  1973  by  R.  H.  Hardin  and  F.  D.  Tappert, 
the  parabolic  equation  (PE)  method  is  an  efficient  means  of  solving  wave  propagation 
problems  [Jensen  et  al.  2000],  The  PE  techniques  have  been  continually  advanced  within 
the  acoustic  modeling  community,  and  they  are  now  the  most  common  method  of 
resolving  range-dependent  underwater  acoustic  propagation  [Jensen  et  al.  2000].  The 
following  PE  derivation  is  closely  aligned  with  that  in  Jensen  et  al  [2000]. 


Beginning  with  the  Helmholtz  equation,  which  describes  acoustic  pressure  as  a 
function  of  range  ( r ),  azimuth  {cp),  and  depth  ( z ),  azimuthal  symmetry  is  assumed  so 
that  the  q>  -coordinate  can  be  ignored.  This  yields  a  description  of  acoustic  pressure  as  a 


function  of  range  and  depth, />(r,  z)  ,  and  the  Helmholtz  equation  becomes: 


d^p  1  dp  d^p 
dr^  r  dr  dz^ 


+  k^n^p  =  0 , 


(4) 


where  =  —  is  the  reference  wavenumber,  n{r,z)  =  — ^ — is  the  index  of  refraction,  and 
c  c(r,z) 

c{r,z)  =  c(z)  -I-  5c{r,z)  .  The  solution  of  the  Helmholtz  equation  is  assumed  to  have  the 
form 

p{r,  z)  =  W (r,  z)H^\k^r) ,  (5) 


where  'P(r,z)  is  the  wave  envelope  function  and  is  a  Hankel  function  of  the  first 
order. 
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By  substituting  equation  (5)  into  (4)  and  making  the  small  angle  approximation, 


— - , 

dr^  “  dr 


(6) 


the  Helmholtz  equation  yields  the  standard  parabolic  equation  as  first  introduced  to 
underwater  acoustics  by  Hardin  and  Tappert: 


2^'^o^  +  ^  +  -^oV-l)^  =  0. 
or  oz 


(7) 


2.  Solutions  to  the  Parabolic  Equation  Using  the  Split  Step  Fourier 
Algorithm 


Split  step  methods  analytically  solve  partial  differential  equations  in 
range-marching  algorithms  with  two  distinct  steps.  In  this  case,  a  solution  is  first 
calculated  in  the  depth  domain  using  a  phase  screen  that  includes  refractive  effects,  and 
then  the  solution  is  advanced  in  the  wavenumber  domain  for  a  homogenous  medium 
including  diffraction.  This  study  uses  the  following  Split  Step  Fourier  algorithm 
described  in  Jensen  et  ah: 


^{r,z)  =  J~  \e 


iAr 

2^ 


JTrri  ^[«^(.o.z)-l]Ar 


'f(r 


O’" 


(8) 


where  is  the  fast  Fourier  transform  (FFT)  from  the  z-domain  to  the  ^z-domain  and 
J"  ■ '  is  the  inverse  transform. 


3.  Set-up  of  the  Computational  Domain 

Three  distinct  boundaries  exist  within  the  computational  domain:  the  ocean 
surface  or  air-ocean  interface,  the  ocean  bottom  or  ocean-seafloor  interface,  and  a 
“sponge”  layer  that  is  used  beneath  the  seafloor  to  prevent  computational  wrap-around  of 
energy  by  the  Fourier  transform  solution. 

Ocean  acoustic  analyses  typically  treat  the  surface  as  a  pressure-release  boundary, 
which  is  represented  by  the  boundary  condition 'F(r,0)  =  0  .  Within  the  PE  model,  this 
boundary  condition  is  achieved  by  creating  an  image  source.  Similar  to  Lloyd’s  mirror, 
this  setup  causes  any  original-source  sound  at  the  surface  to  be  canceled  out  by  the  image 
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source.  By  subtracting  the  image  source  contribution  from  the  original  source 
contribution,  the  surface  boundary  condition  is  met: 

'F(r,0)  =  'F(r,-z,)-'F(r,zJ  =  0,  (9) 

where  z,  is  the  original  source  depth. 


At  the  ocean-seafloor  interface,  the  bottom  sediment  layer  is  important  to 
propagation  because  of  attenuation  and  the  density  change  at  the  interface.  Here,  a 
variable-density  layer  is  included  in  the  numerical  model  to  modify  an  effective  index  of 
refraction,  where  Weff  is  density  dependent.  In  a  real  ocean,  the  density  difference  at  the 
ocean  bottom  is  discontinuous,  but  this  causes  computational  problems  in  solving  the 
numerical  algorithm.  Hence,  we  use 


^  +  Pi)  Pi) 


^z-Z)o^ 


(10) 


V  L  j 

where  D^  is  the  depth  of  the  interface  and  L  is  the  distance  over  which  the  density 
changes  from  to  ,  to  represent  the  density  discontinuity  at  the  interface.  L  is 


defined  such  that  k^L  =  2 . 


Density  of  the  water,  p^,  is  determined  by  the  background  ocean  profile 
(discussed  below),  and  density  of  the  seafloor,  =  1800kg  m'^,  is  representative  of  a 
sandy,  silty  bottom. 


Attenuation  in  the  ocean  bottom  is  treated  by  adding  an  imaginary  part,  ia ,  to  the 
wavenumber, 

CO 


k  =  — \-ia,  a>0 
c 

where  a  is  defined  as  1.11  dB/>..  This  yields  a  complex  index  of  refraction: 

n{x,  zY  -  ^ 


(11) 


rii 

c  \ 

Co 

2 

UoJ 

^c(x,z)J 

(12) 


which  allows  attenuation  to  be  evaluated  in  range  and  depth. 


Below  the  ocean-seafloor  interface,  the  computational  domain  is  treated  with  an 
artificially  attenuating  “sponge”  layer,  in  order  to  absorb  the  periodic  solution  generated 
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by  the  FFT  method.  In  the  periodic  FFT  solution,  energy  at  the  bottom  of  the  solution 
domain  reenters  the  solution  space  at  the  top  of  the  next  cycle.  Artificial  attenuation  in 
the  sponge  layer  of  the  form, 


L{z)  =  exp 


X  exp[- 


(7-7  ^ 
Z  Zjj 

fl 

CCZj,  J 

'j 

(13) 


where  z,,  is  the  depth  at  the  bottom  of  the  computational  field,  negates  this  wrap-around 


effect  by  imposing  a  gradual  loss  of  amplitude  as  z  approaches  the  computational 
bottom  [Colosi  and  Flatte  1996].  Here,  /?  is  the  relative  strength  of  the  loss,  and  a  is 
the  location  relative  to  the  bottom  at  which  loss  is  applied.  As  determined  for  this  study, 
values  of  /3  =  0.04  and  a  =  0.05  stop  acoustic  energy  from  penetrating  beyond  about 
300m  above  the  computational  boundary,  which  is  approximately  1000m  deeper  than  the 
ocean-seafloor  interface. 


4.  Sound  Source 

A  Gaussian  sound  source  is  used  to  closely  align  with  the  far-field,  point-source 
solution  of  the  Helmholtz  equation.  This  source  is  defined  as: 

2  SL/ 

'F(0,z)  =  V^e"^^'''^^10^"”,  (14) 


where  is  the  effective  source  level,  —  is  the  beamwidth,  and  SL  is  the  source  level 

ko 

(195dB).  At  the  surface  boundary,  equation  (9)  is  applied  to  the  Gaussian  sound  source 
to  yield  the  starting  field: 


'F(0,2)  =  V^[e 


(15) 


"2.71  f 

The  unperturbed  wavenumber  can  be  represented  as  = - 2.  terms  of  initial 

Co 

frequency,  f  ^ ,  and  initial  sound  speed,  Cg .  As  an  example,  a  close-up  view  of  a  fiill 
water-column  transmission  loss  diagram  is  shown  in  Figure  5. 
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Figure  5.  Close-up  view  of  transmission  loss  plot  from  250Hz  Gaussian  point  source 
in  an  unperturbed  ocean  at  200m  depth. 


The  250Hz  point  source  in  Figure  5  is  at  200m  and  transmits  sound  energy  in  a 
conical  range  45  degrees  above  and  below  horizontal.  This  45  degree  cone  is  visible 
below  the  horizontal,  with  the  region  void  of  acoustic  energy  represented  by  the  dark  blue 
sliver  on  the  left  hand  side.  In  this  case,  the  sound  source  is  too  close  to  the  surface  for 
the  corresponding  feature  to  be  visible  above  the  source.  The  modeled  transmission  loss 
plot  shows  the  interference  pattern  of  upward  and  downward  acoustic  rays,  as  well  as 
interactions  between  the  bottom  reflections  and  transmitted  energy  at  ranges  greater  than 
10km.  Within  the  first  10km,  the  acoustic  rays  are  steep  enough  that  they  penetrate  into 
the  seafloor,  but  beyond  10km,  energy  is  reflected  off  the  bottom,  causing  the  energy 
visible  above  the  main  ray  loop. 

B.  EDDY  AND  OCEAN  MODELING 
1.  Background  Ocean 

Background  ocean  profiles  for  the  Philippine  Sea  were  derived  from  the  World 
Ocean  Atlas  2001  database.  Annual  average  profiles  for  temperature  and  salinity  from 
the  database  were  used  to  calculate  sound  speed,  potential  density,  buoyancy  frequency, 
and  potential  sound  speed  gradient  profiles.  Figure  6  shows  the  background  ocean 
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profiles  derived  from  the  World  Ocean  Atlas  at  23  °N  125  °E,  which  is  northwest  of  the 
NPAL  tomography  array  (Figure  1)  in  deep  water  greater  than  5000m. 


Figure  6.  The  background  ocean  profiles  used  in  this  study:  (a)  sound  speed, 
(b)  buoyancy  frequency,  (c)  potential  sound  speed  gradient,  (d)  RMS  sound  speed 
variability.  The  profiles  were  derived  from  the  World  Ocean  Atlas  2001  annual 
database  for  23  °N  125  °E.  In  (d),  the  blue  line  is  the  average  of  32  random 
realization  of  the  ocean,  while  the  green  line  is  the  theoretical  profile  for  the 
acoustic  model. 


2.  Eddies 

Eddy-induced  sound  speed  perturbations  are  a  function  of  the  vertical 
displacement,  ,  of  the  water  column  by  the  eddy: 


Sc{r,t)  =  ^{r,t) 


^  dc'^ 
\dzj 


(16) 
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In  order  to  model  the  vertical  displacement  caused  by  the  eddies,  eddy  fields  are 
modeled  as  a  superposition  of  linear  Rossby  waves,  in  accordance  with  Munk  et  al. 
[1995],  This  method  creates  computationally  efficient  eddy  approximations,  which  are 
spatially  and  temporally  correct  in  scale. 


This  method  neglects  the  non-linear  nature  of  ocean  eddies,  but  provides  a 
first-order  approximation  of  the  vertical  displacement  of  the  eddy: 

ax,  y,  Z,  0  =  S  J  ,  (17) 

i=i 

where  aj{k^,  k^)  is  the  Rossby  wave  amplitude,  which  is  a  zero  mean  Gaussian  random 
variable;  IF(z)  is  the  baroclinic  Rossby  mode,  as  depicted  in  Figure  7;  k^,  k^,  and  cr  are 
Rossby  wavenumbers  and  frequency;  and  j  is  the  maximum  Rossby  mode  number 
treated.  In  this  study,  j  =  and  eddies  with  wavelengths  between  25  and  1100km 
and  frequencies  of  0.04  to  3.2  cycles  per  year  were  modeled.  The  spectral  form  used  is: 


SciK) 

Inh 


where  kl  =  k^  +  k^ .  The  normal  mode  equation  is: 


(18) 


dz^  f  ^ 


(19) 


r 


where  =  - 


V 


o-  J 


,  N  is  the  buoyancy  frequency,  and /is  the  Coriolis  parameter. 
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Figure  7.  Normal  mode  eurves. 

Based  on  the  findings  of  Stammer  [1997],  we  assume  isotropie  variability  to 
facilitate  computational  efficiency  in  modeling.  As  shown  in  Figure  4,  two  forms  of  eddy 
wavenumber  spectra  are  used:  that  developed  by  Stammer  [1997]  and  a  Lorentzian  form 
[Munk  et  al.  1995],  both  described  in  Chapter  11. 

As  suggested  by  Wunsch  [1997],  we  treated  80%  of  observed  surface  elevation  as 
due  to  baroclinic  modes,  with  the  baroclinic  energy  distributed  in  the  first  four  modes. 
Therefore,  20%  of  the  total  energy  was  modeled  in  the  barotropic  mode,  which  has  no 
impact  on  sound  speed.  The  remaining  80%  was  modeled  as  baroclinic  energy  and 
further  partitioned  into  Mode  1  (70%),  Mode  2  (15%),  Mode  3  (9%),  and  Mode  4  (6%). 

Based  on  satellite  altimetry,  the  Philippine  Sea  has  approximately  15  cm  root 
mean  square  (RMS)  variability  [Wunsch  1997;  Qui  1999],  as  shown  in  Figures  8  and  9. 
The  Philippine  Sea  variability  is  much  greater  than  most  of  the  eastern  North  Pacific 
Ocean,  where  NPAL  experiments  have  taken  place  in  the  past  and  where  RMS  variability 
is  approximately  5cm.  SSH  variability  of  15cm  RMS  was  used  as  a  normalization  factor, 
and  random  eddy  fields  were  created  to  overlay  on  the  background  ocean.  Examples  of 
random  eddy  fields  are  shown  in  Figure  10,  where  (a)  is  based  on  the  Stammer  spectrum 
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and  appears  to  have  a  correlation  length  around  150km  and  (b)  is  based  on  the  Lorentzian 
spectrum  and  has  a  correlation  length  around  50km,  as  determined  by  the  value  of  k*: 

(20) 


Longitude 


Figures.  Ocean  variability  (y)  from  TOPEX/Poseidon  altimetric  data  from 
December  1992  to  November  1995,  plotted  as  tropical  ocean  (open  circles),  very 
low  energy  areas  with  y  <  6cm  (small  dots),  the  bulk  of  the  oceans  with 
6cm  <y  <15cm  (medium  dots),  and  the  high  energy  areas  with  y  >  15cm.  Large 
bold  circles  are  marked  for  another  purpose  and  can  be  ignored  here.  From 
Stammer  1997. 
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Figure  9.  Map  of  RMS  sea  surface  height  (SSH)  variability  in  the  North  Pacific 
Ocean,  based  on  TOPEX/Poseidon  data  from  October  1992  to  December  1997. 
From  Qui  1999. 


Figure  10.  Map  of  surface  elevation,  modeled  using  the  background  ocean  for 
23°N  125°E  and  random  eddy  fields  from  the  (a)  Stammer  spectrum  and 
(b)  Lorentzian  spectrum. 
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IV.  RESULTS  AND  ANALYSIS 


Eight  different  sets  of  simulations  were  run  with  varying  combinations  of  the 
following  parameters  to  simulate  acoustic  transmission  through  eddy  fields:  acoustic 
frequency  (20  Hz  or  250  Hz),  eddy  wavenumber  spectrum  (Stammer  or  Lorentz),  and 
source  depth  (50m  or  200m),  as  shown  in  Table  1.  The  frequencies  used  were  chosen  to 
represent  possible  tactically-relevant  frequencies;  20Hz  is  near  a  blade-rate  frequency, 
and  250Hz  is  near  machinery  frequencies,  plus  it  is  the  frequency  to  be  used  in  the  NPAL 
experiment.  The  depths  were  chosen  to  simulate  depths  at  which  submarines  might 
transit.  Each  simulation  was  run  for  32  realizations  of  the  ocean,  with  full  water-column 
transmission  loss  (TL)  diagrams  based  on  the  Gaussian  sound  source  (described  in 
Chapter  III)  for  each  realization.  In  addition,  full  water-column  TL  diagrams  were 
generated  for  the  unperturbed  ocean  based  on  the  Gaussian  sound  source. 
Single-receiver,  point-to-point  transmission  loss  curves  from  perturbed  oceans  were 
averaged  and  compared  to  the  unperturbed  ocean  TL  curves  to  assess  the  effects  of  ocean 
mesoscale  features  on  transmission  loss. 

Table  1.  Eight  different  combinations  of  frequency,  spectrum,  and  source  depth 

used  to  generate  TL  curves  through  eddy-field  oceans.  Each  combination  was 
run  32  times. 


Run 

Frequency 

Spectrum 

Source  Depth 

#1 

20  Hz 

Stammer 

50m 

#2 

20  Hz 

Stammer 

200m 

#3 

20  Hz 

Lorentz 

50m 

#4 

20  Hz 

Lorentz 

200m 

#5 

250  Hz 

Stammer 

50m 

#6 

250  Hz 

Stammer 

200m 

#7 

250  Hz 

Lorentz 

50m 

#8 

250  Hz 

Lorentz 

200m 
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A. 


BACKGROUND  OCEAN  AND  EDDY  FIELD  TRANSMISSION  LOSS 


Each  of  the  eight  simulations  generated  one  full  water-column  TL  diagram  for  the 
unperturbed  ocean  and  32  plots  for  the  eddy  realizations.  Examples  of  each  combination 
are  shown  in  Figures  11-18.  The  TL  diagrams  on  the  left  are  examples  of  transmission 
loss  through  an  unperturbed  ocean,  while  the  TL  diagrams  on  the  right  are  examples  of 
transmission  loss  through  a  modeled  eddy-field.  Cold  colors  indicate  areas  of  high 
transmission  loss,  while  warm  colors  indicate  lower  transmission  loss.  These  TL 
diagrams  do  not  include  cylindrical  spreading,  and  the  seafloor  is  depicted  by  the  red  line 
at  5250m. 

Convergence  zones  are  defined  as  regions  in  the  deep  ocean  where  sound  rays, 
refracted  from  the  depths,  return  to  the  surface.  The  rays  are  focused  at  or  near  the 
surface  in  successive  intervals  [Commander,  Navy  Warfare  Development  2006].  If  these 
intervals  are  viewed  in  plan  view,  they  appear  as  concentric  circles  around  a  source, 
which  are  referred  to  as  convergence  zone  (CZ)  rings.  In  the  vertical  cross-sections 
shown  in  Figures  11-18,  convergence  zones  appear  as  the  brighter  loops  of  green,  yellow 
and  orange  through  the  entire  water  column.  CZs  have  horizontal  lengths  on  order  50km, 
with  focusing  regions  (rings)  near  the  surface.  For  instance,  in  Figure  11a,  the 
convergence  zone  rings  are  approximately  at  60km,  120km,  180km,  250km,  310km,  and 
370km. 

Figures  11-18  demonstrate  some  of  the  common  physical  processes  caused  by  the 
eddy  fields.  In  all  the  figures  (11-18),  there  is  noticeable  vertical  and  horizontal  shifting 
of  the  CZs  and  spreading  of  CZs.  Specifically,  Figure  18b  is  an  example  of  considerable 
shifting  of  the  CZ,  as  acoustic  rays  get  refracted  at  the  third  CZ  ring,  near  175km  at  500m 
depth,  creating  a  secondary  set  of  CZ  loops  and  partially  filling  in  the  shadow  zones  (dark 
blue  areas)  that  existed  in  the  unperturbed  case  (Figure  18a).  To  clarify,  the  CZ  rings  in 
18a  are  near  65km,  130km,  190km,  250km,  310km,  and  375km.  In  Figure  I8b,  the  first 
two  CZ  rings  are  similar:  60km,  120km,  but  at  the  third  ring  (~  175km),  acoustic  rays 
from  the  near-source  side  of  the  transmission  “beam”  are  diverted,  creating  two  beams  of 
transmission.  The  original  beam  has  CZ  rings  at  250km,  310km,  and  380km,  while  the 

newly  created  beam  has  CZ  rings  at  210km,  260km,  325km,  and  385km. 
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In  a  few  cases,  the  eddy  fields  alter  TL  by  driving  energy  into  the  seafloor  (as 
shown  in  15b)  and  even  into  the  surface  duct  (as  in  Figures  13b,  14b,  and  17b).  Energy 
driven  into  the  seafloor  causes  significantly  damped  transmission  after  the  interaction, 
indicated  by  darker  blue  colors  in  the  TL  plot.  For  example,  in  Figure  15b,  the  first 
bottom  interaction  is  around  35km,  after  which  there  is  less  turquoise  and  more  of  the 
darker  blues,  indicating  increased  transmission  loss  as  compared  to  the  unperturbed 
simulation  (Figure  15a).  Energy  trapped  in  surface  ducts  is  depicted  by  a  yellow  or  green 
horizontal  line  of  energy  near  the  surface,  for  example  in  Figure  13b  from  the  origin  to 
150km.  In  Figure  1 4b,  this  occurs  at  the  fourth  and  fifth  CZ  rings,  about  250km  and 
325km;  trapping  is  evident  in  Figure  I7b  near  250km.  Scintillation  is  quantified  by  the 
varying  intensities  between  realizations,  so  its  effects  are  not  depicted  in  these  figures. 
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Figure  11.  (a)  Full  water-column  transmission  loss  diagram  for  unperturbed  ocean: 

frequency  20  Hz,  source  depth  50m.  (b)  As  in  (a),  but  with  ocean  perturbed  by 
Stammer  eddy  field.  In  (b),  the  eddy  field  causes  spreading  and  shifting  of  the 
CZs. 


Figure  12.  (a)  Full  water-column  transmission  loss  diagram  for  unperturbed  ocean: 

frequency  20  Hz,  source  depth  200m.  (b)  As  in  (a),  but  with  ocean  perturbed  by 
Stammer  eddy  field.  In  (b),  the  eddy  field  causes  spreading  and  shifting  of  the 
CZs. 


24 


Figure  13.  (a)  Full  water-column  transmission  loss  diagram  for  unperturbed  ocean: 

frequency  250  Hz,  source  depth  50m.  (b)  As  in  (a),  but  with  ocean  perturbed  by 
Stammer  eddy  field.  In  (b),  the  eddy  field  causes  energy  to  become  trapped  in  the 
surface  duct  near  the  origin.  There  is  also  spreading  and  shifting  of  the  CZs 
created  by  energy  that  escapes  the  surface  duct. 


(a)  Range  (km)  (b)  Range  (km) 

Figure  14.  (a)  Full  water-column  transmission  loss  diagram  for  unperturbed  ocean: 

frequency  250  Hz,  source  depth  200m.  (b)  As  in  (a),  but  with  ocean  perturbed  by 
Stammer  eddy  field.  In  (b),  the  eddy  field  causes  spreading  and  shifting  of  the 
CZs,  as  well  as  some  trapped  energy  in  the  surface  duct  near  250  and  325  km. 
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Figure  15.  (a)  Full  water-column  transmission  loss  diagram  for  unperturbed  ocean: 

frequency  20  Hz,  source  depth  50m.  (b)  As  in  (a),  but  with  ocean  perturbed  by 
Lorentz  eddy  field.  In  (b),  the  eddy  field  causes  significant  energy  to  be  trapped 
in  the  seafloor,  as  shown  by  the  darker  blue  colors  in  (b),  especially  after  the  first 
bottom  interaction  around  35km.  The  eddy  field  also  causes  spreading  and 
shifting  of  the  CZs. 


Figure  16.  (a)  Full  water-column  transmission  loss  diagram  for  unperturbed  ocean: 

frequency  20  Hz,  source  depth  200m.  (b)  As  in  (a),  but  with  ocean  perturbed  by 
Lorentz  eddy  field.  In  (b),  the  eddy  field  causes  spreading  and  shifting  of  the 
CZs. 


26 


Figure  17.  (a)  Full  water-column  transmission  loss  diagram  for  unperturbed  ocean: 

frequency  250  Hz,  source  depth  50m.  (b)  As  in  (a),  but  with  ocean  perturbed  by 
Lorentz  eddy  field.  In  (b),  the  eddy  field  causes  energy  to  be  trapped  in  the 
surface  duct  around  250  km,  as  well  as  spreading  and  shifting  of  the  CZs. 


Figure  18.  (a)  Full  water-column  transmission  loss  diagram  for  unperturbed  ocean: 

frequency  250  Hz,  source  depth  200m.  (b)  As  in  (a),  but  with  ocean  perturbed  by 
Lorentz  eddy  field.  In  (b),  the  eddy  field  causes  significant  spreading  and  shifting 
of  the  CZ  energy,  which  partially  fills  the  shadow  zones.  The  major  deviation 
between  the  unperturbed  and  the  eddy  realization  occurs  around  175km  at  about 
500m  depth  in  (b). 
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A  quick  visual  examination  of  the  figures  above  reveals  differenees  in  the  CZ 
struetures  for  high  and  low  frequency  sources.  Focusing  on  the  unperturbed  plots,  two 
distinct  processes  are  causing  the  different  struetures  visible  in  Figures  1  la- 18a.  The  first 
is  evaneseent  waves,  whieh  cause  striations  of  TL  in  the  shadow  zones.  For  the  20Hz 
source,  the  evanescent  waves  deeay  more  slowly,  leaving  more  energy  (brighter  eolors) 
in  the  shadow  zones.  For  the  250Hz  source,  these  evaneseent  waves  deeay  mueh  more 
rapidly,  so  the  shadow  zones  are  predominantly  free  of  acoustie  energy  [Flatte  1979]. 
Secondly,  the  lower  frequency  acoustic  energy  is  more  affected  by  diffraction,  leading  to 
mueh  more  diffuse  and  broadened  eonvergenee  zone  loops  [Flatte  1979].  This  is 
espeeially  evident  in  comparing  Figures  11a  and  13a.  The  eonvergenee  zone  loops  for 
the  higher  frequency  (Figure  13a)  are  mueh  more  distinct  because  they  are  less  affeeted 
by  diffraetion  and  the  associated  evaneseent  waves  deeay  more  quickly.  Also,  the  higher 
frequeney  source  is  better  foeused  near  the  surface,  creating  more  distinet  CZ  rings 
[Flatte  1979]. 

Additionally,  the  lower  frequency  has  a  longer  wavelength,  so  it  has  a  larger-seale 
interference  pattern,  whieh  leads  to  fewer  interferenee  bands  or  lobes  in  the  eonvergenee 
zone  loops  [Flatte  1979].  This  is  espeeially  evident  with  the  deeper  souree  depth. 
Comparing  Figures  12a  and  14a,  the  interference  lobes  for  the  lower  frequeney  souree  are 
not  deteetable,  while  in  Figure  14a,  at  least  nine  lobes  are  visually  identifiable. 

B.  EFFECTS  ON  TRANSMISSION  LOSS 

For  eaeh  of  the  eight  model  runs,  point-to-point  transmission  loss  curves  to  a 
single  reeeiver  at  50m  depth  were  calculated  for  the  unperturbed  ocean  (<l)o(r))  and  for 

each  realization  of  the  eddy  fields  ( 0,.(r) ,  where  i  is  the  realization  number).  The  TL 

eurves  were  smoothed  with  range  using  a  Butterworth  filter  in  order  to  eliminate  the  high 
variability  associated  with  multipath  interference  patterns.  The  smoothing  serves  to 
eliminate  the  small-seale  shifts  from  the  interference  pattern  and  to  focus  statistics  on  the 
large-scale  ehanges  in  intensity.  The  smoothing  process  is  equivalent  to  averaging  over 
range,  whieh  is  algebraieally  similar  to  averaging  over  frequency,  so  smoothing  these 
curves  makes  them  equivalent  to  TL  eurves  one  might  see  with  a  broadband  signal 
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[Harrison  and  Harrison  1995].  Hence,  even  though  single  frequencies  are  modeled,  the 
results  presented  here  are  reasonable  estimations  what  might  occur  during  the  NPAL 
experiment  where  broadband  signals  will  be  used. 
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Figure  19.  An  example  of  unsmoothed  XL  curve  (blue)  showing  much  variability  due 
to  multipath  interference.  The  red  curve  is  smoothed  using  the  Butterworth  filter. 


After  smoothing,  the  32  curves  from  each  run  were  averaged  with  each  other  to 
create  a  composite  eddy  XL  curve  for  each  scenario: 

(  1  i=32  A 


O£(r)  =  101ogj„| 


(21) 


V  i  J 

where  is  the  composite  eddy-induced  transmission  loss  in  dB  and  Ii(r)  is  the  acoustic 
intensity  at  range  r  for  each  realization  i.  These  composite  curves  were  then  differenced 
with  the  unperturbed  ocean  curve 

0,(r)  =  0^(r)-101ogl0(/o(r)),  (22) 

where  0^(r)  is  the  bias  curve  and  Io(r)  is  the  unperturbed  acoustic  intensity  at  range  r. 


The  Oj(r)  curves  are  presented  below  in  Figures  20-32  as  the  top  panel  in  (a)  and  (b) 

and  represent  the  expected  increase  or  decrease  of  acoustic  intensity  through  an  eddy 
field. 
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Root  mean  square  variability  for  the  32  perturbed  curves  was  calculated  using  the 
composite  eddy  TL  curve  as  the  mean: 


f  1  ‘=^2 


f  1  /=32  ^2 


(23) 


where  =  101ogj(,/(r,/) .  represents  the  possible  variability  of  sound 

intensity  around  the  mean  TL  curve;  it  is  depicted  in  the  bottom  panel  of  Figures  20-32  in 
(a)  and  (b). 


In  order  to  draw  valid  summary  statistics  from  these  two  curves 
(<I)j(r)  and  ),  a  threshold  value  above  a  reasonable  ambient  noise  level  was  set 

for  each  unperturbed  TL  curve.  This  created  an  index  array  of  all  range  positions  that 
were  above  the  threshold  value.  The  and  values  associated  with  this  index 

were  used  to  calculate  the  summary  statistics.  This  method  emphasizes  the  regions  where 
acoustic  energy  is  significant,  which  are  also  regions  of  high  variability,  and  eliminates 
TL  bias  and  variability  where  the  sound  level  would  be  too  low  to  be  detectable. 

The  summary  statistics  used  to  analyze  the  trends  among  the  three  variables  are 
significant  bias  ( ),  which  is  a  quantification  of  the  shift  in  detectable  sound  caused  by 

the  eddy  field;  and  significant  RMS  which  is  a  measure  of  the  variability  of 

sound  through  the  eddy  field.  takes  into  account  three  of  the  mechanisms  by  which 

eddies  alter  TL  -  shifting  CZs,  driving  energy  into  the  seafloor,  and  trapping  energy  in  a 
surface  duct.  incorporates  the  fourth  mechanism  -  scintillation.  Summary 

statistics  are  presented  in  Table  2,  and  they  will  be  discussed  further  in  the  paragraphs  1-3 
below. 


Expected  uncertainties  for  individual  estimates  of  Oj(r)  and  are  of  the 

order 

statistics  reduces  this  error  by  about  a  factor  2.  Therefore,  the  uncertainty  or  error 
associated  with  the  summary  statistics  is  about  8%.  This  becomes  important  when 
discussing  the  significance  of  and  values. 


or  about  17%;  however,  range  averaging  to  obtain  the  significant  summary 
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Table  2.  Summary  statistics:  <I)^  is  a  measure  of  the  expected  decibel  shift  of 

acoustic  energy  through  an  eddy  ocean  versus  the  unperturbed  case  (measured 
only  where  the  signal  has  high  intensity  above  the  given  threshold).  is  the 

square-root  of  significant  variance,  a  measure  of  the  variance  of  the  32 
realizations  around  their  mean  (calculated  for  above-threshold  signal). 


Run 

Freq  (Hz), 
Spectra, 

Source  Depth  (m) 

O,  (dB) 

#1 

20,  Stammer,  50 

0.59 

2.17 

#2 

20,  Stammer,  200 

0.80 

3.33 

#3 

20,  Lorentz,  50 

0.59 

2.36 

#4 

20,  Lorentz,  200 

0.88 

3.02 

#5 

250,  Stammer,  50 

5.33 

6.18 

#6 

250,  Stammer,  200 

5.69 

6.58 

#7 

250,  Lorentz,  50 

4.37 

5.50 

#8 

250,  Lorentz,  200 

5.46 

5.66 

1.  Frequency  Dependency 

For  each  frequency,  four  combinations  of  spectra  and  source  depth  were 
computed,  as  shown  in  Table  1,  which  allows  for  four  comparisons  of  frequency 
tendencies  on  transmission  loss.  Figures  20  through  23  show  these  comparisons,  where 
the  top  two  panels  are  0^(r)  ,  and  the  bottom  two  panels  show  ^^^(r) .  The  plots  on  the 

left  are  for  the  20Hz  source  (a),  and  the  plots  on  the  right  are  for  the  250Hz  source  (b). 
Summary  statistics  are  shown  in  Table  2,  and  Table  3  is  a  comparison  of  the  statistics  for 
the  two  different  frequencies  modeled. 
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Figure  20.  (a)  0^(r)  and  ,  using  Stammer  speetrum,  frequeney  20  Hz,  souree 

depth  50m,  reeeiver  depth  50m.  Statisties  for  run  #1  show  0^=0. 59dB  and 
<I)^5=2.17dB.  (b)  As  in  (a),  but  with  frequeney  250  Hz.  Statistics  for  run  #5 
show  <l)^=5.33dB  and  <l)^^5=6.18dB. 


Figure  21.  (a)  <I>j(r)  and  ,  using  Stammer  spectrum,  frequency  20  Hz,  source 

depth  200m,  receiver  depth  50m.  Statistics  for  run  #2  show  (I>^=0.80dB  and 
0^^=3.33dB.  (b)  As  in  (a),  but  with  frequency  250  Hz.  Statistics  for  run  #6 
show  <l)^=5.69dB  and  <l)^^5=6.58dB. 
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Figure  22.  (a)  Oj(r)  and  ,  using  Lorentzian  spectrum,  frequency  20  Hz, 

source  depth  50m,  receiver  depth  50m.  Statistics  for  run  #3  show  =0.59dB 
and  =2.36dB.  (b)  As  in  (a),  but  with  frequency  250  Hz.  Statistics  for  run  #7 
show  <I)j=4.37dB  and  <I)j^5=5.50dB. 


Figure  23.  (a)  and  ,  using  Lorentzian  spectrum,  frequency  20  Hz, 

source  depth  200m,  receiver  depth  50m.  Statistics  for  run  #4  show  =0.88dB 
and  =3.02dB.  (b)  As  in  (a),  but  with  frequency  250  Hz.  Statistics  for  run  #8 
show  <I)j=5.46dB  and  <I)^5=5.66dB. 
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Table  3.  Comparison  of  statistics  for  low  frequency  (20Hz)  transmissions  versus 

higher  frequency  (250Hz)  transmissions  through  the  modeled  eddy  fields. 
Consistently,  the  acoustic  energy  from  the  higher  frequency  source  experienced 
more  and  than  the  low  frequency  acoustic  energy. 


Run 

Freq  (Hz),  Spectra, 
Source  Depth  (m) 

d)ft  (dB) 

(dB) 

Run 

Freq  (Hz),  Spectra, 
Source  Depth  (m) 

d)ft  (dB) 

(dB) 

#1 

20,  Stammer,  50 

0.59 

2.17 

#5 

250,  Stammer,  50 

5.33 

6.18 

#2 

20,  Stammer,  200 

0.80 

3.33 

#6 

250,  Stammer,  200 

5.69 

6.58 

#3 

20,  Lorentz,  50 

0.59 

2.36 

#7 

250,  Lorentz,  50 

4.37 

5.50 

#4 

20,  Lorentz,  200 

0.88 

3.02 

#8 

250,  Lorentz,  200 

5.46 

5.66 

Figures  20-23  and  Table  3  show  that  d)^  values  for  20Hz  are  less  than  one 

decibel,  while  the  corresponding  values  for  250Hz  are  about  five  decibels.  This  result  is 
somewhat  expected  because  lower  frequencies  are  more  stable,  and  the  acoustic  energy  is 
not  as  affected  by  mesoscale  features  as  the  higher  frequency  energy.  For  the  higher 
frequency,  there  is  more  energy  refracted  into  the  shadow  zones,  which  allows  detectable 
acoustic  energy  where  there  previously  had  been  none.  The  lower  frequency  energy  also 
was  not  as  susceptible  to  refraction  and  trapping  in  the  surface  duct  because  the  lower 
acoustic  frequencies  have  longer  wavelengths  which  are  not  trapped  as  easily  in  the 
surface  layer. 

From  Figures  20-23  and  Table  3,  the  d)^^  values  for  20Hz  are  two  to  three 

decibels,  and  the  corresponding  values  for  250Hz  are  five  to  six  decibels.  This  result  was 
also  expected  due  to  greater  stability  of  the  low  frequency  signal.  The  variability 
measurement  is  a  quantification  of  scintillation  occurring,  which  is  caused  by  multipath 
interference.  The  high  frequency  signal  is  more  likely  to  break  into  micro-rays,  which 
are  the  primary  source  of  multipath  interference. 

Tactically,  this  dependence  on  frequency  indicates  that  less  variability  is  expected 
with  low  frequency  sources.  In  other  words,  acoustic  prediction  for  low  frequency  tones 
does  not  require  mesoscale  features  for  accurate  ranges. 
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2.  Dependency  on  Spectrum 

For  each  spectrum,  four  combinations  of  frequency  and  source  depth  were 
computed,  as  shown  in  Table  1,  which  allows  for  four  comparisons  of  spectral  effects  on 
transmission  loss.  Figures  24  through  27  show  these  comparisons,  with  the  top  two 
panels  being  Oj(r) ,  and  the  bottom  two  panels  showing  .  The  plots  on  the  left 

are  for  the  Stammer  eddy  spectrum  (a),  and  the  plots  on  the  right  are  for  the  Lorentzian 
eddy  spectrum  (b).  The  Stammer  spectrum  has  correlation  lengths  around  150km,  while 
the  Lorentzian  spectrum  has  correlation  lengths  around  50km.  Summary  statistics  are 
shown  in  Table  2,  and  a  side-by-side  comparison  of  the  statistics  for  the  two  spectral 
cases  is  shown  in  Table  4. 


_  1 

0  50  100  150  200  2S0  300  350  400 

Range  (km) 

■  l.il  ll 

II: 

■  ■  ■  V 

0  50  100  150  200  250  300  350400 

Range  <ltfn) 


Range  (km> 


Figure  24.  (a)  <I>j(r)  and  ,  using  Stammer  spectrum,  frequency  20  Hz,  source 

depth  50m,  receiver  depth  50m.  Statistics  for  run  #1  show  0^=0. 59dB  and 
0^^=2.17dB.  (b)  As  in  (a),  but  with  Lorentzian  spectrum.  Statistics  for  run  #3 
show  <I)j=0.59dB  and 
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Figure  25.  (a)  <I>^(r)  and  ,  using  Stammer  spectrum,  frequency  20  Hz,  source 

depth  200m,  receiver  depth  50m.  Statistics  for  run  #2  show  0^=0.80dB  and 
®flMs^3.33dB.  (b)  As  in  (a),  but  with  Lorentzian  spectrum.  Statistics  for  run  #4 
show  0^=0. 88dB  and  <I)j,jj^5=3.02dB. 


^3^  R»»9*(km)  Range(km) 

Figure  26.  (a)  <I>j(r)  and  ,  using  Stammer  spectrum,  frequency  250  Hz, 

source  depth  50m,  receiver  depth  50m.  Statistics  for  run  #5  show  <I)j=5.33dB 
and  =6. 18dB.  (b)  As  in  (a),  but  with  Lorentzian  spectrum.  Statistics  for  run 
#7  show  0^=4.37dB  and  Og^5=5.50dB. 
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Figure  27.  (a)  <I>j(r)  and  ,  using  Stammer  spectrum,  frequency  250  Hz, 

source  depth  200m,  receiver  depth  50m.  Statistics  for  run  #6  show  <I)j=5.69dB 
and  =6.58dB.  (b)  As  in  (a),  but  with  Lorentzian  spectrum.  Statistics  for  run 
#8  show  (l>^=5.46dB  and  (l>g^5=5.66dB. 

Table  4.  Comparison  of  statistics  for  Stammer  spectrum  versus  Lorentzian 

spectrum  for  the  modeled  eddy  fields.  The  two  spectra  produce  similar  numbers 
for  Oj  and  in  all  combinations  of  variables. 


Run 

Freq  (Hz),  Spectra, 
Source  Depth  (m) 

(dB) 

(dB) 

Run 

Freq  (Hz),  Spectra, 
Source  Depth  (m) 

0)^  (dB) 

(dB) 

#1 

20,  Stammer,  50 

0.59 

2.17 

#3 

20,  Lorentz,  50 

0.59 

2.36 

#2 

20,  Stammer,  200 

0.80 

3.33 

#4 

20,  Lorentz,  200 

0.88 

3.02 

#5 

250,  Stammer,  50 

5.33 

6.18 

#7 

250,  Lorentz,  50 

4.37 

5.50 

#6 

250,  Stammer,  200 

5.69 

6.58 

#8 

250,  Lorentz,  200 

5.46 

5.66 

Figures  24-27  and  Table  4  show  that  values  for  both  Stammer  and  Lorentzian 

eddies  are  less  than  one  decibel  for  the  low  frequency  case,  with  the  difference  between 
corresponding  cases  for  the  two  spectra  being  only  a  few  tenths  of  a  decibel.  For  the 
higher  frequency  examples,  the  two  spectra  cause  similar  values  for  ,  right  around  five 

decibels  for  both  spectra.  Based  on  expected  uncertainties  for  ,  a  few  tenths  of  a 

decibel  is  not  a  significant  difference.  An  insignificant  trend  between  the  two  spectra  is  a 

surprising  result  considering  how  different  the  sea  surface  height  realizations  looked 

(Figure  10).  The  correlation  length  for  the  Lorentzian  eddies  is  approximately  50km,  and 
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the  correlation  length  for  the  Stammer  eddies  appears  to  be  approximately  150km,  yet 
these  seem  to  have  little  to  do  with  the  changes  in  acoustic  transmission. 

From  Figures  24-27  and  Table  4,  the  values  for  both  Stammer  and 

Lorentzian  eddies  are  two  to  three  decibels  at  the  low  frequency  and  five  to  six  decibels 
at  the  higher  frequency.  Again,  based  on  expected  uncertainties  for  ,  these  results 

are  not  a  significant  difference.  The  results  suggest  that  the  eddy  wavenumber  spectrum 
has  little  effect  on  the  acoustic  transmissions.  Comuelle  and  Howe’s  research  [1987] 
revealed  that  loop  resonance  structure  and  scales  are  important  for  travel  time  variance, 
yet  this  study  suggests  that  intensity  and  arrival  time  are  sensitive  to  different  scales. 
Until  the  issue  of  intensity  sensitivity  is  better  resolved,  tactical  decisions  about  model 
resolution  cannot  accurately  be  made. 

3.  Source  Depth  Dependency 

For  each  source  depth,  four  combinations  of  spectra  and  frequency  were 
computed,  as  shown  in  Table  1,  which  allows  for  four  comparisons  of  source  depth 
tendencies  on  transmission  loss.  Figures  28  through  3 1  show  these  comparisons,  with  the 
top  two  panels  being  ,  and  the  bottom  two  panels  showing  .  The  plots  on  the  left 

are  for  the  source  at  50m  (a),  and  the  plots  on  the  right  are  for  the  source  at  200m  (b). 
Summary  statistics  are  shown  in  Table  2,  and  Table  5  is  a  comparison  of  the  statistics  for 
the  two  source  depths. 
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Figure  28.  (a)  0^(r)  and  ,  using  Stammer  spectrum,  frequency  20  Hz,  source 

depth  50m,  receiver  depth  50m.  Statistics  for  run  #1  show  0^=0.5906  and 
<I)^5=2.17dB.  (b)  As  in  (a),  but  with  source  depth  200m.  Statistics  for  run  #2 
show  <l)^=0.80dB  and  <l)^^5=3.33dB. 


Figure  29.  (a)  <I>^(r)  and  ,  using  Stammer  spectrum,  frequency  250  Hz, 

source  depth  50m,  receiver  depth  50m.  Statistics  for  run  #5  show  0^=5.33dB 
and  <I)j^5=6.18dB.  (b)  As  in  (a),  but  with  source  depth  200m.  Statistics  for  run 
#6  show  <I)^=5.69dB  and  <I)^j=6.58dB. 
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Figure  30.  (a)  Oj(r)  and  ,  using  Lorentzian  spectrum,  frequency  20  Hz, 

source  depth  50m,  receiver  depth  50m.  Statistics  for  run  #3  show  =0.59dB 
and  <l)^5=2.36dB.  (b)  As  in  (a),  but  with  source  depth  200m.  Statistics  for  run 
#4  show  (I>^=0.88dB  and  (I>^5=3.02dB. 


Figure  31.  (a)  0^(r)  and  ,  using  Lorentzian  spectrum,  frequency  250  Hz, 

source  depth  50m,  receiver  depth  50m.  Statistics  for  run  #7  show  =4.37dB 
and  <l)^5=5.50dB.  (b)  As  in  (a),  but  with  source  depth  200m.  Statistics  for  run 
#8  show  0^=5.46dB  and  Og^5=5.66dB. 
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Table  5.  Comparison  of  statistics  for  shallow  source  depth  (50m)  versus  deeper 

source  depth  (200m)  through  the  modeled  eddy  fields.  The  depths  produce 
similar  orders  of  magnitude  for  <I)^  and  in  all  combinations  of  variables,  but 

the  deeper  source  depth  consistently  has  slightly  higher  values  for  both  and 


Run 

Freq  (Hz),  Spectra, 
Source  Depth  (m) 

(dB) 

(dB) 

Run 

Freq  (Hz),  Spectra, 
Source  Depth  (m) 

(dB) 

(dB) 

#1 

20,  Stammer,  50 

0.59 

2.17 

#2 

20,  Stammer,  200 

0.80 

3.33 

#3 

20,  Lorentz,  50 

0.59 

2.36 

#4 

20,  Lorentz,  200 

0.88 

3.02 

#5 

250,  Stammer,  50 

5.33 

6.18 

#6 

250,  Stammer,  200 

5.69 

6.58 

#7 

250,  Lorentz,  50 

4.37 

5.50 

#8 

250,  Lorentz,  200 

5.46 

5.66 

Figures  28-31  and  Table  5  show  that  values  for  the  50m  source  are  about  six 

tenths  of  a  decibel  for  the  low  frequency  case  and  are  between  four  and  five  decibels  for 
the  high  frequency  case.  With  the  200m  source,  the  values  for  the  low  frequency 

case  are  just  under  one  decibel  with  the  high  frequency  case  experiencing  greater  than 
five  decibels  .  The  source  at  200m  consistently  has  higher  values  for  ,  regardless 

of  frequency.  Again,  the  expected  uncertainty  in  this  calculation  is  about  8%,  so  these 
differences  in  are  more  significant  than  the  uncertainty.  Two  physical  processes  are 

influencing  this  bias:  bottom  interaction  is  causing  negative  bias,  and  CZ  broadening  is 
shifting  energy  into  the  shadow  zones,  which  is  a  positive  bias.  The  shallow  source  has 
more  transmission  loss  to  attenuation  in  the  seafloor,  and  the  negative  bias  values 
associated  with  that  TL  keep  the  shallow-source  near  zero,  despite  the  positive 

influence  of  CZ  broadening.  For  the  deeper  source,  the  bottom  loss  is  not  as  significant, 
so  the  CZ  broadening  is  the  dominant  process  influencing  the  deep-source  . 

From  Figures  28-31  and  Table  5,  the  values  for  the  50m  source  are  about 

two  decibels  for  the  low  frequency  case  and  about  between  five  and  six  decibels  for  the 
high  frequency  case.  The  deeper  source  depth  has  values  greater  than  three 

decibels  for  the  low  frequency  case  and  near  six  decibels  for  the  higher  frequency  case. 

Despite  the  frequency  dependence,  the  deeper  source  depth  consistently  experienced 
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higher  variability  than  the  shallower  source  depth.  The  differences  in  are 

significant  because  the  expected  uncertainty  is  approximately  8%,  or  less  than  one 
decibel. 

To  sum  up  this  section,  factors  that  increase  both  and  are  higher 

frequency  and  deeper  source  depth.  The  spectrum  used  for  the  simulations  does  not 
appear  to  have  a  significant  effect  on  acoustic  variability. 


C.  OUT-OF-PLANE  SCATTERING 


There  was  not  sufficient  time  in  this  study  to  simulate  the  range-dependent 
horizontal  deflections  that  occur  as  sound  transmits  through  eddies;  however,  the  effect 
on  beam-forming  and  target  angle  localization  is  significant,  so  some  theory  is  discussed 
here.  Much  of  this  section  is  modified  from  Flatte  and  Colosi  [2008],  but  it  is  thoroughly 
discussed  in  other  literature  as  well  [Munk  1980,  Weinberg  and  Clark  1980,  and  Munk  et 
al.  1995]. 


Using  a  ray-based  physical  model  to  simulate  acoustic  transmission,  sound  speed 
is  represented  as 


c(x,  y,  z)  =  c(z)  +  Sc(x,  y,  z) , 


(24) 


where  Jc  is  a  small  perturbation  when  compared  to  c(z) ,  the  mean  sound  speed.  The 
fractional  sound  speed  change  is  defined  as 


^0 


(25) 


Horizontal  angular  deviations  of  rays  due  to  refraction  through  sound  speed 
fluctuations  can  be  derived  from  the  ray  equations: 


do  ^  d/u  _  , 

dx  dy 


(26) 


where  -/j.'  scales  as  —  and  is  the  random  gradient  of  fractional  sound  speed  (in  the 
y-direction).  The  deviation  after  traveling  a  distance  dx  is  approximately  dO  --ju'dx . 
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In  order  to  convert  this  angular  deviation  into  horizontal  out-of-plane  distance,  the 
integration  must  first  be  calculated  in  time,  where  the  angle  of  the  ray  at  range  R  is: 

0(R)  ju\x)dx.  (27) 

Here,  the  integration  in  time  over  the  perturbed  ray  path  s(t)  is  approximated  as  the 
integration  over  range  of  the  unperturbed  ray  path.  This  approximation  leads  to  a  useful 
form  of  variance  in  ray  angle: 


This  form  of  the  integration  uses  Lh  and  Ly  to  represent  the  correlation  lengths  along  the 
ray,  ,  and  assumes  that  the  correlation  length  of  the  sound  speed  gradient  is 

negligible  compared  to  range,  R.  As  shown  above,  the  RMS  angle  variance  in  the 
horizontal  grows  as  the  square-root  of  range.  Taking  typical  eddy-scale  values, 
=  2x10“^  and  Ly  =  50  km ,  the  RMS  horizontal  angle  is  approximately  equal  to: 

-  0.8°  (R/500  km)^  (29) 

at  a  range  of  500km. 

When  compared  to  an  internal  wave  horizontal  angle  deviation  representative  of 
similar  ocean  conditions  ((//2)  =4.0x10'^  and  =  10  km ), 

{d^f^  ~  0.25°  (R/500  km)^ ,  (30) 

the  angular  deviation  caused  by  eddies  is  significantly  greater,  nearly  by  a  factor  of  three. 


Using  the  small  angle  approximation  to  calculate  range  error  from  angular 
deviation,  is  approximately  equal  to  ,  as  shown  in  Figure  32,  where 


S0=(e:f 


Combining  the  formulas  above. 


(  R  ^2 

y  Lfj  )  v500  km 


(31) 
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and  using  typical  eddy-scale  values,  the  RMS  range  error  at  500  km  due  to  ocean  eddies 
is 


dy  =  l  km 


R 


(32) 


500  km  j 

When  compared  to  an  internal  wave  range  error  representative  of  similar  ocean 
conditions  as  above, 


5y  =  2.3  km 


500  km 


(33) 


the  range  error  at  500  km  caused  by  eddies  is  significantly  greater,  by  a  factor  of  three. 


R 


Figure  32.  Angular  deviation,  <56’ ,  where  R  is  unperturbed  range  and  Syis  range 
error  due  to  angular  deviation. 


The  scientific  community  has  put  much  effort  into  quantifying  internal  wave 
effects  on  acoustic  transmission.  The  above  demonstration  shows  that  ocean  mesoscale 
features,  such  as  eddies,  are  at  least  as  important  as  internal  waves  in  causing  horizontal 
variations  in  acoustic  transmission. 
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V.  CONCLUSION 


A.  CURRENT  EXPERIMENT 

A  two-dimensional  parabolic  equation  model  was  built  to  simulate  acoustic 
transmission  through  mesoscale  environments  in  the  Philippine  Sea.  The  model  used 
three  essential  parts.  The  background  ocean  was  based  on  climatological  ocean  profiles 
for  the  Philippine  Sea.  The  eddy  fields  were  constructed  using  two  wavenumber  spectra: 
Stammer  spectrum  with  correlation  length  150km  and  Lorentzian  spectrum  with 
correlation  length  50km.  These  eddy  fields  were  then  overlaid  on  the  background  ocean. 
Next,  the  acoustic  portion  of  the  model  was  incorporated  using  parabolic-equation 
propagation  of  a  Gaussian  sound  source.  Global  altimetry  data  revealed  that 
approximately  15cm  of  RMS  variability  in  sea  surface  elevation  exists  in  the  Philippine 
Sea.  This  sea  surface  variability  represents  local  oceanic  kinetic  energy.  Based  on 
Wunsch’s  vertical  partitioning  of  oceanic  energy,  20%  of  kinetic  energy  was  modeled  in 
the  barotropic  modes  and  80%  in  the  baroclinic  modes.  The  baroclinic  energy  was  then 
further  partitioned  into  the  first  four  baroclinic  modes.  Eddies  were  simulated  as  a  linear 
superposition  of  Rossby  waves  of  the  first  four  baroclinic  modes.  This  model  was  used 
to  analyze  acoustic  propagation  fluctuations  through  the  background  (climatological) 
ocean  and  the  perturbed  (eddy-field)  ocean. 

Eight  different  combinations  of  acoustic  frequency  (20Hz  and  250Hz),  eddy 
wavenumber  spectra  (Stammer  or  Eorentzian),  and  source  depth  (50m  or  200m)  were 
simulated  through  this  model.  For  each  combination,  the  model  was  run  for  32  random 
realizations  of  eddy  fields.  For  each  of  these  realizations,  a  point-to-point  transmission 
loss  (TE)  curve  was  calculated  from  source  depth  (50m  or  200m)  to  a  receiver  at  50m. 
This  TL  curve  was  also  calculated  for  the  unperturbed  ocean.  These  curves  were 
smoothed  in  range  using  a  Butterworth  filter  to  eliminate  the  variability  associated  with 
multipath  interference.  From  the  32  perturbed  curves,  a  composite  eddy-field  TL  curve 
was  derived,  and  then  this  curve  was  differenced  from  the  unperturbed  curve  to  assess  the 
acoustic  variability  caused  by  the  mesoscale  ocean  features. 
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Several  physical  phenomena  influence  the  patterns  of  transmission  loss  that  are 
examined  in  this  study.  Evanescent  waves  occur  outside  of  the  convergence  zones,  and 
they  decay  much  less  with  low  frequency  signals.  This  causes  the  cluttered  appearance 
of  the  low  frequency  full  water-column  TL  plots.  Low  frequency  signals  also  experience 
more  diffraction  than  higher  frequency  transmissions,  leading  to  broader  convergence 
zones.  Conversely,  high  frequency  signals  experience  stronger  focusing  near  the  surface, 
so  these  simulations  show  much  stronger  CZ  rings  near  the  surface.  The  high  frequency 
signals  have  a  smaller-scale  interference  pattern  due  to  their  shorter  wavelengths.  This 
creates  the  appearance  of  several  lobes  of  energy  within  the  main  CZ  beam  in  the  TL 
plots. 

Four  physical  mechanisms  caused  by  eddies  influence  TL.  is  a  measure  of 

the  convergence  zone  broadening,  driving  energy  into  the  seafloor,  and  trapping  energy 
in  the  surface  duct.  <I)^j  measures  the  fourth  mechanism  -  scintillation. 

The  computer  simulations  revealed  higher  frequency  sound  transmitting  through 
the  eddy  fields  and  deeper  source  depths  cause  the  largest  amount  of  acoustic  variability, 
both  in  the  bias  (O^)  and  in  the  RMS  variability  (O^^).  These  results  suggest  that 

acoustic  variability  through  eddies  is  strongly  dependent  on  frequency  and  weakly 
dependent  on  source  depth.  Somewhat  surprisingly,  the  wavenumber  spectrum  used  to 
simulate  eddies  had  no  significant  effect  on  the  acoustic  variability  beyond  that  which 
was  already  caused  by  the  frequency  or  source  depth. 

The  frequency  dependence  is  so  strong  because  the  convergence  zone  shifting  is 
highly  frequency  dependent.  Low  frequencies  already  have  broad  CZs  due  to  diffraction, 
so  the  eddy  fields  do  not  have  as  much  an  impact,  while  for  the  higher  frequency  signals, 
the  CZ  are  significantly  broadened  due  to  the  eddies.  In  addition,  the  lower  frequency 
signal  was  less  likely  to  get  trapped  in  the  surface  duct,  which  was  another  factor  that 
limited  its  .  The  higher  frequency  is  more  likely  to  split  into  micro-rays,  which 

increases  its  . 
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The  spectrum  dependence  is  insignificant,  but  revealed  an  area  for  further  study. 
Clearly,  travel  time  and  intensity  are  more  sensitive  to  different  spectra  of  ocean  features, 
so  more  clarity  is  required  before  any  modeling  decisions  are  based  on  eddy  wavenumber 
spectra. 

Source  depth  weakly  impacts  the  trends  for  and  .  The  CZ  spreading 

causes  a  positive  bias,  while  energy  trapped  into  the  seafloor  causes  a  negative  bias. 
Increased  scintillation  and  surface  duct  trapping  do  not  appear  to  influence  d)^  or  d)^^ 

trends  as  they  relate  to  source  depth. 

The  simulations  for  this  study  revealed  fluctuations  of  sound  intensity  of  5-20dB, 
which  are  significant  to  naval  applications  for  the  impact  to  detection  methods  and 
capabilities.  For  the  upcoming  NPAL  field  experiment  in  the  Philippine  Sea,  the 
significance  applies  more  to  the  usefulness  of  the  tomographic  data  collected  and  how  to 
best  process  that  data  at  the  conclusion  of  the  experiment’s  field  work.  It  also  validates 
the  plan  to  map  the  mesoscale  features  of  the  Philippine  Sea,  because  a  significant 
advantage  can  be  gained  through  accurate  knowledge  of  the  mesoscale. 

A  final  area  of  consideration  is  the  horizontal  out-of-plane  refraction  caused  by 
eddies.  This  can  cause  a  horizontal  shift  on  order  one  degree  at  500km,  which  is  about 
7km  of  localization  error  at  500km. 

B.  FUTURE  RESEARCH 

Eddy  effects  on  horizontal  out-of-plane  scattering  are  important  because  they 
affect  beam-forming  and  target  localization  accuracy.  Additional  computer  simulation  of 
these  effects  is  necessary  to  improving  acoustic  operations  in  a  mesoscale-rich 
environment. 

Some  questions  that  this  study  uncovered  also  warrant  further  research.  Despite 
the  large  difference  in  sea  surface  elevation  correlation  lengths  (and  corresponding  eddy 
size)  caused  by  the  two  different  spectra  (Stammer  and  Lorentz),  there  was  little 
discemable  difference  in  the  acoustic  variability  caused  by  each  spectrum.  Further 
investigation  into  the  variables  that  impact  acoustic  variability  should  be  pursued. 
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The  model  developed  here  should  be  compared  to  existing  acoustic  models  to 
determine  the  impact  of  the  mesoscale  environment  beyond  what  is  already  being 
modeled.  Additionally,  these  computer  simulations  need  to  be  compared  to  actual 
acoustic  data  to  see  if  nature  reinforces  the  principles  suggested  by  the  simulations. 
However,  the  comparison  becomes  difficult  because  it  is  nearly  impossible  to  separate 
eddy  effects  from  other  mesoscale  effects  such  as  internal  waves  in  real  ocean  data. 

A  fourth  area  of  research  is  in  practical  application  of  these  findings.  Tactically, 
how  does  the  crew  of  a  submarine  exploit  the  ocean  variability  to  make  their  own  ship 
less  detectable  while  simultaneously  expanding  their  ability  to  detect  enemy  vessels? 

A  useful  outcome  of  the  NPAL  experiment  will  be  actual  measurements  of  the 
proportional  mode  energy  and  baroclinic  structures  of  eddies,  as  well  as  the  important 
scales  of  eddies  for  acoustics.  This  is  important  to  compare  to  the  theory  presented  by 
Cornuelle  and  Howe  [1987]  that  waves  with  lengths  comparable  to  the  ray  loop  length 
are  most  acoustically  significant. 


The  horizontal  out-of-plane  scattering  has  important  impacts  on  beam-forming 
capabilities.  The  horizontal  coherence  length  is  dictated  by  and  the  local 

wavenumber: 


k  2 

/tfl 


(34) 


where  (6^)  is  the  RMS  angle  deviation  (described  in  Chapter  IV,  Section  C)  and  is 

the  local  wavenumber.  Based  on  the  angle  deviation  presented  here,  what  are  the  effects 
to  beam-forming  coherence? 


Preliminary  results  presented  here  demonstrate  that  ocean  eddies,  especially  the 
powerful  ones  in  the  Philippine  Sea,  can  cause  significant  fluctuations  to  acoustic 
transmission  loss  in  the  ocean.  With  further  research,  acoustic  prediction  models  can 
evolve  from  range-independent  to  high-resolution,  eddy-resolving,  range-dependent 
models,  which  will  be  important  in  truly  understanding  the  ocean  battlespace  and 
exploiting  that  knowledge  to  naval  tactical  advantage. 
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